Mapping canopy traits over Québec using airborne and spaceborne imaging spectroscopy.

The advent of new spaceborne imaging spectrometers offers new opportunities for ecologists to map vegetation traits at global scales. However, to date most imaging spectroscopy studies exploiting satellite spectrometers have been constrained to the landscape scale. In this paper we present a new method to map vegetation traits at the landscape scale and upscale trait maps to the continental level, using historical spaceborne imaging spectroscopy (Hyperion) to derive estimates of leaf mass per area, nitrogen, and carbon concentrations of forests in Québec, Canada. We compare estimates for each species with reference field values and obtain good agreement both at the landscape and continental scales, with patterns consistent with the leaf economic spectrum. By exploiting the Hyperion satellite archive to map these traits and successfully upscale the estimates to the continental scale, we demonstrate the great potential of recent and upcoming spaceborne spectrometers to benefit plant biodiversity monitoring and conservation efforts.

Seasonal trends in leaf-level photosynthetic capacity and water use efficiency in a North American Eastern deciduous forest and their impact on canopy-scale gas exchange.

Vegetative transpiration (E) and photosynthetic carbon assimilation (A) are known to be seasonally dynamic, with changes in their ratio determining the marginal water use efficiency (WUE). Despite an understanding that stomata play a mechanistic role in regulating WUE, it is still unclear how stomatal and nonstomatal processes influence change in WUE over the course of the growing season. As a result, limited understanding of the primary physiological drivers of seasonal dynamics of canopy WUE remains one of the largest uncertainties in earth system model projections of carbon and water exchange in temperate deciduous forest ecosystems. We investigated seasonal patterns in leaf-level physiological, hydraulic, and anatomical properties, including the seasonal progress of the stomatal slope parameter (g1 ; inversely proportional to WUE) and the maximum carboxylation rate (Vcmax ). Vcmax and g1 were seasonally variable; however, their patterns were not temporally synchronized. g1 generally showed an increasing trend until late in the season, while Vcmax peaked during the midsummer months. Seasonal progression of Vcmax was primarily driven by changes in leaf structural, and anatomical characteristics, while seasonal changes in g1 were most strongly related to changes in Vcmax and leaf hydraulics. Using a seasonally variable Vcmax and g1 to parameterize a canopy-scale gas exchange model increased seasonally aggregated A and E by 3% and 16%, respectively.

The effect of the vertical gradients of photosynthetic parameters on the CO2 assimilation and transpiration of a Panamanian tropical forest.

Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies. We compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama and then quantified the impact of the observed gradients on simulated canopy-scale CO2 and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy-scale simulations of water vapor and CO2 exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top-of-canopy, leaf-level water-use efficiency was markedly higher at the top-of-canopy, and the decrease in maximum carboxylation rate from the top-of-canopy to the ground was less than TBM assumptions. The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within-individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species-rich, complex tropical forests.

Integrating plant physiology into simulation of fire behavior and effects.

Wildfires are a global crisis, but current fire models fail to capture vegetation response to changing climate. With drought and elevated temperature increasing the importance of vegetation dynamics to fire behavior, and the advent of next generation models capable of capturing increasingly complex physical processes, we provide a renewed focus on representation of woody vegetation in fire models. Currently, the most advanced representations of fire behavior and biophysical fire effects are found in distinct classes of fine-scale models and do not capture variation in live fuel (i.e. living plant) properties. We demonstrate that plant water and carbon dynamics, which influence combustion and heat transfer into the plant and often dictate plant survival, provide the mechanistic linkage between fire behavior and effects. Our conceptual framework linking remotely sensed estimates of plant water and carbon to fine-scale models of fire behavior and effects could be a critical first step toward improving the fidelity of the coarse scale models that are now relied upon for global fire forecasting. This process-based approach will be essential to capturing the influence of physiological responses to drought and warming on live fuel conditions, strengthening the science needed to guide fire managers in an uncertain future.

Short-term variation in leaf-level water use efficiency in a tropical forest.

The representation of stomatal regulation of transpiration and CO2 assimilation is key to forecasting terrestrial ecosystem responses to global change. Given its importance in determining the relationship between forest productivity and climate, accurate and mechanistic model representation of the relationship between stomatal conductance (gs ) and assimilation is crucial. We assess possible physiological and mechanistic controls on the estimation of the g1 (stomatal slope, inversely proportional to water use efficiency) and g0 (stomatal intercept) parameters, using diurnal gas exchange surveys and leaf-level response curves of six tropical broadleaf evergreen tree species. g1 estimated from ex situ response curves averaged 50% less than g1 estimated from survey data. While g0 and g1 varied between leaves of different phenological stages, the trend was not consistent among species. We identified a diurnal trend associated with g1 and g0 that significantly improved model projections of diurnal trends in transpiration. The accuracy of modeled gs can be improved by accounting for variation in stomatal behavior across diurnal periods, and between measurement approaches, rather than focusing on phenological variation in stomatal behavior. Additional investigation into the primary mechanisms responsible for diurnal variation in g1 will be required to account for this phenomenon in land-surface models.

Leaves as bottlenecks: The contribution of tree leaves to hydraulic resistance within the soil-plant-atmosphere continuum.

Within vascular plants, the partitioning of hydraulic resistance along the soil-to-leaf continuum affects transpiration and its response to environmental conditions. In trees, the fractional contribution of leaf hydraulic resistance (Rleaf ) to total soil-to-leaf hydraulic resistance (Rtotal ), or fRleaf (=Rleaf /Rtotal ), is thought to be large, but this has not been tested comprehensively. We compiled a multibiome data set of fRleaf using new and previously published measurements of pressure differences within trees in situ. Across 80 samples, fRleaf averaged 0.51 (95% confidence interval [CI] = 0.46-0.57) and it declined with tree height. We also used the allometric relationship between field-based measurements of soil-to-leaf hydraulic conductance and laboratory-based measurements of leaf hydraulic conductance to compute the average fRleaf for 19 tree samples, which was 0.40 (95% CI = 0.29-0.56). The in situ technique produces a more accurate descriptor of fRleaf because it accounts for dynamic leaf hydraulic conductance. Both approaches demonstrate the outsized role of leaves in controlling tree hydrodynamics. A larger fRleaf may help stems from loss of hydraulic conductance. Thus, the decline in fRleaf with tree height would contribute to greater drought vulnerability in taller trees and potentially to their observed disproportionate drought mortality.

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean (Glycine max) nested association mapping population.

Photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500-2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R2 from 0.56 to 0.96 and root mean square error approximately <10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis.

Late-day measurement of excised branches results in uncertainty in the estimation of two stomatal parameters derived from response curves in Populus deltoides Bartr. × Populus nigra L.

Many terrestrial biosphere models depend on an understanding of the relationship between stomatal conductance and photosynthesis. However, unlike the measurement of photosynthetic parameters, such as the maximum carboxylation capacity, where standard methods (e.g., CO2 response or ACi curves) are widely accepted, a consensus method for empirically measuring parameters representing stomatal response has not yet emerged. Most models of stomatal response to environment represent stomatal conductance as being bounded by a lower intercept parameter (g0), and linearly scaled based on a multivariate term described by the stomatal slope parameter (g1). Here we employ the widely used Unified Stomatal Optimization model, to test whether g1 and g0 parameters are impacted by the choice of measurement method, either on an intact branch or a cut branch segment stored in water. We measured paired stomatal response curves on intact and excised branches of a hybrid poplar clone (Populus deltoides Bartr. × Populus nigra L. OP367), measured twice over a diurnal period. We found that predawn branch excision did not significantly affect measured g0 and g1 when measured within 4 h of excision. Measurement in the afternoon resulted in significantly higher values of g1 and lower values of g0, with values changing by 55% and 56%, respectively. Excision combined with afternoon measurement resulted in a marked effect on parameter estimates, with g1 increasing 89% from morning to afternoon and a 25% lower g1 for cut branches than those measured in situ. We also show that in hybrid poplar the differences in parameter estimates obtained from plants measured under different conditions can directly impact models of canopy function, reducing modeled transpiration by 18% over a simulated 12.5-h period. Although these results are only for a single isohydric woody species, our findings suggest that stomatal optimality parameters may not remain constant throughout the day.

An improved representation of the relationship between photosynthesis and stomatal conductance leads to more stable estimation of conductance parameters and improves the goodness-of-fit across diverse data sets.

Stomata play a central role in surface-atmosphere exchange by controlling the flux of water and CO2 between the leaf and the atmosphere. Representation of stomatal conductance (gsw ) is therefore an essential component of models that seek to simulate water and CO2 exchange in plants and ecosystems. For given environmental conditions at the leaf surface (CO2 concentration and vapor pressure deficit or relative humidity), models typically assume a linear relationship between gsw and photosynthetic CO2 assimilation (A). However, measurement of leaf-level gsw response curves to changes in A are rare, particularly in the tropics, resulting in only limited data to evaluate this key assumption. Here, we measured the response of gsw and A to irradiance in six tropical species at different leaf phenological stages. We showed that the relationship between gsw and A was not linear, challenging the key assumption upon which optimality theory is based-that the marginal cost of water gain is constant. Our data showed that increasing A resulted in a small increase in gsw at low irradiance, but a much larger increase at high irradiance. We reformulated the popular Unified Stomatal Optimization (USO) model to account for this phenomenon and to enable consistent estimation of the key conductance parameters g0 and g1 . Our modification of the USO model improved the goodness-of-fit and reduced bias, enabling robust estimation of conductance parameters at any irradiance. In addition, our modification revealed previously undetectable relationships between the stomatal slope parameter g1 and other leaf traits. We also observed nonlinear behavior between A and gsw in independent data sets that included data collected from attached and detached leaves, and from plants grown at elevated CO2 concentration. We propose that this empirical modification of the USO model can improve the measurement of gsw parameters and the estimation of plant and ecosystem-scale water and CO2  fluxes.

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities.

As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiagvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

Reducing model uncertainty of climate change impacts on high latitude carbon assimilation.

The Arctic-Boreal Region (ABR) has a large impact on global vegetation-atmosphere interactions and is experiencing markedly greater warming than the rest of the planet, a trend that is projected to continue with anticipated future emissions of CO2 . The ABR is a significant source of uncertainty in estimates of carbon uptake in terrestrial biosphere models such that reducing this uncertainty is critical for more accurately estimating global carbon cycling and understanding the response of the region to global change. Process representation and parameterization associated with gross primary productivity (GPP) drives a large amount of this model uncertainty, particularly within the next 50 years, where the response of existing vegetation to climate change will dominate estimates of GPP for the region. Here we review our current understanding and model representation of GPP in northern latitudes, focusing on vegetation composition, phenology, and physiology, and consider how climate change alters these three components. We highlight challenges in the ABR for predicting GPP, but also focus on the unique opportunities for advancing knowledge and model representation, particularly through the combination of remote sensing and traditional boots-on-the-ground science.

Rapid estimation of photosynthetic leaf traits of tropical plants in diverse environmental conditions using reflectance spectroscopy.

Tropical forests are one of the main carbon sinks on Earth, but the magnitude of CO2 absorbed by tropical vegetation remains uncertain. Terrestrial biosphere models (TBMs) are commonly used to estimate the CO2 absorbed by forests, but their performance is highly sensitive to the parameterization of processes that control leaf-level CO2 exchange. Direct measurements of leaf respiratory and photosynthetic traits that determine vegetation CO2 fluxes are critical, but traditional approaches are time-consuming. Reflectance spectroscopy can be a viable alternative for the estimation of these traits and, because data collection is markedly quicker than traditional gas exchange, the approach can enable the rapid assembly of large datasets. However, the application of spectroscopy to estimate photosynthetic traits across a wide range of tropical species, leaf ages and light environments has not been extensively studied. Here, we used leaf reflectance spectroscopy together with partial least-squares regression (PLSR) modeling to estimate leaf respiration (Rdark25), the maximum rate of carboxylation by the enzyme Rubisco (Vcmax25), the maximum rate of electron transport (Jmax25), and the triose phosphate utilization rate (Tp25), all normalized to 25°C. We collected data from three tropical forest sites and included leaves from fifty-three species sampled at different leaf phenological stages and different leaf light environments. Our resulting spectra-trait models validated on randomly sampled data showed good predictive performance for Vcmax25, Jmax25, Tp25 and Rdark25 (RMSE of 13, 20, 1.5 and 0.3 µmol m-2 s-1, and R2 of 0.74, 0.73, 0.64 and 0.58, respectively). The models showed similar performance when applied to leaves of species not included in the training dataset, illustrating that the approach is robust for capturing the main axes of trait variation in tropical species. We discuss the utility of the spectra-trait and traditional gas exchange approaches for enhancing tropical plant trait studies and improving the parameterization of TBMs.

Detection of the metabolic response to drought stress using hyperspectral reflectance.

Drought is the most important limitation on crop yield. Understanding and detecting drought stress in crops is vital for improving water use efficiency through effective breeding and management. Leaf reflectance spectroscopy offers a rapid, non-destructive alternative to traditional techniques for measuring plant traits involved in a drought response. We measured drought stress in six glasshouse-grown agronomic species using physiological, biochemical, and spectral data. In contrast to physiological traits, leaf metabolite concentrations revealed drought stress before it was visible to the naked eye. We used full-spectrum leaf reflectance data to predict metabolite concentrations using partial least-squares regression, with validation R2 values of 0.49-0.87. We show for the first time that spectroscopy may be used for the quantitative estimation of proline and abscisic acid, demonstrating the first use of hyperspectral data to detect a phytohormone. We used linear discriminant analysis and partial least squares discriminant analysis to differentiate between watered plants and those subjected to drought based on measured traits (accuracy: 71%) and raw spectral data (66%). Finally, we validated our glasshouse-developed models in an independent field trial. We demonstrate that spectroscopy can detect drought stress via underlying biochemical changes, before visual differences occur, representing a powerful advance for measuring limitations on yield.

Spectroscopy outperforms leaf trait relationships for predicting photosynthetic capacity across different forest types.

Leaf trait relationships are widely used to predict ecosystem function in terrestrial biosphere models (TBMs), in which leaf maximum carboxylation capacity (Vc,max ), an important trait for modelling photosynthesis, can be inferred from other easier-to-measure traits. However, whether trait-Vc,max relationships are robust across different forest types remains unclear. Here we used measurements of leaf traits, including one morphological trait (leaf mass per area), three biochemical traits (leaf water content, area-based leaf nitrogen content, and leaf chlorophyll content), one physiological trait (Vc,max ), as well as leaf reflectance spectra, and explored their relationships within and across three contrasting forest types in China. We found weak and forest type-specific relationships between Vc,max and the four morphological and biochemical traits (R2 ≤ 0.15), indicated by significantly changing slopes and intercepts across forest types. By contrast, reflectance spectroscopy effectively collapsed the differences in the trait-Vc,max relationships across three forest biomes into a single robust model for Vc,max (R2 = 0.77), and also accurately estimated the four traits (R2 = 0.75-0.94). These findings challenge the traditional use of the empirical trait-Vc,max relationships in TBMs for estimating terrestrial plant photosynthesis, but also highlight spectroscopy as an efficient alternative for characterising Vc,max and multitrait variability, with critical insights into ecosystem modelling and functional trait ecology.

A best-practice guide to predicting plant traits from leaf-level hyperspectral data using partial least squares regression.

Partial least squares regression (PLSR) modelling is a statistical technique for correlating datasets, and involves the fitting of a linear regression between two matrices. One application of PLSR enables leaf traits to be estimated from hyperspectral optical reflectance data, facilitating rapid, high-throughput, non-destructive plant phenotyping. This technique is of interest and importance in a wide range of contexts including crop breeding and ecosystem monitoring. The lack of a consensus in the literature on how to perform PLSR means that interpreting model results can be challenging, applying existing models to novel datasets can be impossible, and unknown or undisclosed assumptions can lead to incorrect or spurious predictions. We address this lack of consensus by proposing best practices for using PLSR to predict plant traits from leaf-level hyperspectral data, including a discussion of when PLSR is applicable, and recommendations for data collection. We provide a tutorial to demonstrate how to develop a PLSR model, in the form of an R script accompanying this manuscript. This practical guide will assist all those interpreting and using PLSR models to predict leaf traits from spectral data, and advocates for a unified approach to using PLSR for predicting traits from spectra in the plant sciences.

Source:sink imbalance detected with leaf- and canopy-level spectroscopy in a field-grown crop.

The finely tuned balance between sources and sinks determines plant resource partitioning and regulates growth and development. Understanding and measuring metabolic indicators of source or sink limitation forms a vital part of global efforts to increase crop yield for future food security. We measured metabolic profiles of Cucurbita pepo (zucchini) grown in the field under carbon sink limitation and control conditions. We demonstrate that these profiles can be measured non-destructively using hyperspectral reflectance at both leaf and canopy scales. Total non-structural carbohydrates (TNC) increased 82% in sink-limited plants; leaf mass per unit area (LMA) increased 38% and free amino acids increased 22%. Partial least-squares regression (PLSR) models link these measured functional traits with reflectance data, enabling high-throughput estimation of traits comprising the sink limitation response. Leaf- and canopy-scale models for TNC had R2 values of 0.93 and 0.64 and %RMSE of 13 and 38%, respectively. For LMA, R2 values were 0.91 and 0.60 and %RMSE 7 and 14%; for free amino acids, R2 was 0.53 and 0.21 with %RMSE 20 and 26%. Remote sensing can enable accurate, rapid detection of sink limitation in the field at the leaf and canopy scale, greatly expanding our ability to understand and measure metabolic responses to stress.

Seasonal trends in photosynthesis and leaf traits in scarlet oak.

Understanding seasonal variation in photosynthesis is important for understanding and modeling plant productivity. Here, we used shotgun sampling to examine physiological, structural and spectral leaf traits of upper canopy, sun-exposed leaves in Quercus coccinea Münchh (scarlet oak) across the growing season in order to understand seasonal trends, explore the mechanisms underpinning physiological change and investigate the impact of extrapolating measurements from a single date to the whole season. We tested the hypothesis that photosynthetic rates and capacities would peak at the summer solstice, i.e., at the time of peak photoperiod. Contrary to expectations, our results reveal a late-season peak in both photosynthetic capacity and rate before the expected sharp decrease at the start of senescence. This late-season maximum occurred after the higher summer temperatures and vapor pressure deficit and was correlated with the recovery of leaf water content and increased stomatal conductance. We modeled photosynthesis at the top of the canopy and found that the simulated results closely tracked the maximum carboxylation capacity of Rubisco. For both photosynthetic capacity and modeled top-of-canopy photosynthesis, the maximum value was therefore not observed at the summer solstice. Rather, in each case, the measurements at and around the solstice were close to the overall seasonal mean, with values later in the season leading to deviations from the mean by up to 41 and 52%, respectively. Overall, we found that the expected Gaussian pattern of photosynthesis was not observed. We conclude that an understanding of species- and environment-specific changes in photosynthesis across the season is essential for correct estimation of seasonal photosynthetic capacity.

Beyond ecosystem modeling: A roadmap to community cyberinfrastructure for ecological data-model integration.

In an era of rapid global change, our ability to understand and predict Earth's natural systems is lagging behind our ability to monitor and measure changes in the biosphere. Bottlenecks to informing models with observations have reduced our capacity to fully exploit the growing volume and variety of available data. Here, we take a critical look at the information infrastructure that connects ecosystem modeling and measurement efforts, and propose a roadmap to community cyberinfrastructure development that can reduce the divisions between empirical research and modeling and accelerate the pace of discovery. A new era of data-model integration requires investment in accessible, scalable, and transparent tools that integrate the expertise of the whole community, including both modelers and empiricists. This roadmap focuses on five key opportunities for community tools: the underlying foundations of community cyberinfrastructure; data ingest; calibration of models to data; model-data benchmarking; and data assimilation and ecological forecasting. This community-driven approach is a key to meeting the pressing needs of science and society in the 21st century.

Multi-hypothesis comparison of Farquhar and Collatz photosynthesis models reveals the unexpected influence of empirical assumptions at leaf and global scales.

Mechanistic photosynthesis models are at the heart of terrestrial biosphere models (TBMs) simulating the daily, monthly, annual and decadal rhythms of carbon assimilation (A). These models are founded on robust mathematical hypotheses that describe how A responds to changes in light and atmospheric CO2 concentration. Two predominant photosynthesis models are in common usage: Farquhar (FvCB) and Collatz (CBGB). However, a detailed quantitative comparison of these two models has never been undertaken. In this study, we unify the FvCB and CBGB models to a common parameter set and use novel multi-hypothesis methods (that account for both hypothesis and parameter variability) for process-level sensitivity analysis. These models represent three key biological processes: carboxylation, electron transport, triose phosphate use (TPU) and an additional model process: limiting-rate selection. Each of the four processes comprises 1-3 alternative hypotheses giving 12 possible individual models with a total of 14 parameters. To broaden inference, TBM simulations were run and novel, high-resolution photosynthesis measurements were made. We show that parameters associated with carboxylation are the most influential parameters but also reveal the surprising and marked dominance of the limiting-rate selection process (accounting for 57% of the variation in A vs. 22% for carboxylation). The limiting-rate selection assumption proposed by CBGB smooths the transition between limiting rates and always reduces A below the minimum of all potentially limiting rates, by up to 25%, effectively imposing a fourth limitation on A. Evaluation of the CBGB smoothing function in three TBMs demonstrated a reduction in global A by 4%-10%, equivalent to 50%-160% of current annual fossil fuel emissions. This analysis reveals a surprising and previously unquantified influence of a process that has been integral to many TBMs for decades, highlighting the value of multi-hypothesis methods.